Supersolidity, entropy and frustration 
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We study the properties of t — t' — V model of hard-core bosons on the triangular lattice that 
can be realized in optical lattices. By mapping to the spin-1/2 XXZ model in a field, we determine 
the phase diagram of the t — V model where the supersolid characterized by the ordering pattern 
(x,x, —2a;') ( "ferrimagnetic" or SS A) is a ground state for chemical potential /j, > 3V . By turning 
on either temperature or t' at half-filling (fj, = 3V), we find a first order transition from SS A to the 
elusive supersolid characterized by the (x, — x, 0) ordering pattern ( "antiferromagnetic" or SS C). In 
addition, we find a large region where a superfluid phase becomes a solid upon raising temperature 
at fixed chemical potential. This is an analog of the Pomeranchuk effect driven by the large entropic 
effects associated with geometric frustration on the triangular lattice. 
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PACS numbers: 75.10. Jm,05. 30. Jp, 67.40.Kh,74.25.Dw 

Supersolidity is one of the most intriguing properties 
of matter. In that state, matter can flow without viscos- 
ity, like in a superfluid, yet atoms arc located at regular 
positions: Translation and U (1) symmetry are broken 
simultaneously. It was originally proposed [l[ that this 
state could exist in 4 He. While such a supersolid state 
may have been observed, Q it is likely that the relevant 
mechanism for 4 He is disorder not zero point vacan- 
cies as first envisioned. 

To observe supersolidity without disorder, one can load 
ultracold bosonic atoms into optical lattices [J. In- 
deed, Bose-Einstein condensation (BEC) of chromium 
atoms in an optical trapping potential |5[ has already 
been observed, making it likely that supersolid phases 
on such lattices can eventually be achieved. Tempera- 
ture is clearly an extremely relevant parameter for these 
experiments. Q 

One of the most promising lattices to observe super- 
solid phases is the triangular lattice where supersolidity 
appears as a result of geometric frustration, from a kind 
of order- by-disorder mechanism @, 0, H, & El], U | • Su- 
persolidity in other two-dimensional lattice models has 
been predicted theoretically, but the triangular lattice 
offers a particularly rich and interesting phase diagram 
in a lattice that is simple to realize. For example, it has 
been proposed [T3| that second-neighbor hopping may 
induce the intriguing particle-hole symmetric supersolid 
C phase, (so-called "antiferromagnetic" supersolid). It 
has been conjectured fl2| that the transition between su- 
persolid C and other phases, such as supersolid A ("fer- 
rimagnetic" supersolid), could occur through a critical 
point with emergent degrees of freedom that cannot be 
described by the standard Landau theory. 

In this paper, we obtain detailed phase diagrams show- 
ing that a particle-hole symmetric supersolid phase C can 
indeed be stabilized by both next-nearest-neighbor hop- 
ping and by finite temperature effects. In addition, the 
frustration associated with the triangular lattice ampli- 



fies entropic effects, leading to a wide range of parameters 
where one can observe supcrfluid-solid-liquid transitions 
as temperature is raised at constant chemical potential. 
On the square lattice [3] > this sequence of transitions oc- 
curs in an extremely narrow range of chemical potentials. 
This phenomenon is an analog of the Pomeranchuk effect 
in 3 He where liquid (not superfluid)-solid-liquid transi- 
tions are observed by increasing T at fixed pressure. 

Model: We consider hard core bosons (infinite on- 
site repulsion) on a triangular lattice, with both near- 
est neighbor (nn) hopping and repulsion (i, V) and next 
nearest neighbor (mm) hopping (t ) 



H = — tijo\aj + h.c + V riiTij — ji > 



(i) 



<ij> 



where each lattice site can be occupied by or 1 boson 
(rii = 0, 1) rii = a\a,i, and fi is the chemical potential. In 
the above restricted Hilbert space, the model (p} can be 
mapped to the S =1/2 XXZ model in a field (h) 

n = vY, sts z 3 - E ^ s i s t + h - c - - h E s i > ( 2 ) 

<ij> i,j i 

where h = [i — 3V. In this language, supersolid 
(SS) ordering corresponds to spins having their x — y 
component aligned ferromagnetically (superfluid (SF)) 
along with their z-component also ordered but at 
non zero wave vector inside the first Brillouin zone 
(solid (S)). A phase without ordering but non-zero 
z-componcnt and zero x — y component corresponds 
to the normal fluid (NF). Fully polarized up (down) 
spins corresponds to Full (Empty) lattice. The order 
parameter for the solid (staggered magnetization in 
spin language, staggered density in boson language) 
is defined with the help of the three sublattice mag- 
netizations {Sf = TH-1/2J = 1,2,3) as (IBj M s = 
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measures the solid order, i.e. a periodicity longer than 
that of the underlying lattice. 
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Method: The Berezinskii-Kosterlitz-Thouless (BKT) 
transitions and the SS-C phase (in the t/V — ► limit) 
are normally out of reach for simple mean-field theo- 
ries. Using large enough clusters however, Self-consistent 
Cluster Mean Field Theory (SCMFT) can overcome some 
of these deficiencies. We argue that, while not perfectly 
accurate, SCMFT [lj§] is an extremely efficient way of 
exploring vast uncharted territory in the phase diagram. 
More refined studies can then improve the accuracy of 
phase boundaries in a second stage. We briefly describe 
the method and then demonstrate its accuracy by com- 
paring with known results. 

A cluster '1' with a finite number of sites (shaded re- 
gion of the inset of Fig. 1) is embedded in the effective 
field of its surroundings. In other words we consider the 
following cluster C spin Hamiltonian TL S 

H s = Y,H ij +Y,KSt+h.c.+Y,htSt-hY,Si (3) 

i,jeC ieC ieC ieC 

where h~ and hf are the effective fields of the surround- 
ings. TL S needs to be diagonalized with the following 
self-consistency conditions: 

v /+ = E ^ < s 7 /+ >> h ' = E ^ < s j > ( 4 ) 

j 3 

where j indicates neighbor of site i and prime over E indi- 
cates that sites j inside the cluster are excluded. Average 
values in Eq.((l|) are obtained from H s . 

Validity of the approach: To assess the accuracy of 
SCMFT, we first show that it reproduces quite accurately 
the phase diagram obtained from the most reliable ap- 
proaches. From now on, we discuss the results mostly in 
the bosonic language. For the t — V model, consider a 
cluster 'V of three sites shown as a shaded area in the 
inset of Fig. 1. We measure t,p,, and temperature T in 
units of V, defining t = t/V,fi = fi/V, and f = T/V. 

We display the zero-temperature phase diagram in 
Figffl This phase diagram is very close to the phase dia- 
gram obtained by Quantum Monte Carlo (QMC) meth- 
ods in Refs. @, Q. In the simplest mean-field approach, 
[13] the supersolid region at Jl = 3 is much too large, ex- 
tending to (t) MF = 0.5 compared with (t) QMC = 0.124 
in Ref. 0. Here, we obtain 0.216, closer to QMC. Also, 
in our approach, the maximum extent of the solid re- 
gion, t = 0.22, is overestimated by only 10% compared 
with the QMC result 0.195. In Fig. 1 we show the spin- 
odals Jls (t) and JIsf (t) between which metastable phases 
or coexistence of SF and S may occur. 

In supersolid A (SS A, Jl > 3), the density on three con- 
secutive sites follows the "ferrimagnetic" ordering pat- 
terns < rii — j >=< S" >= (.x, x, —2x'). In supersolid B 
(SS B, Jl < 3) the pattern is (— x, — x, 2a;') with x ^ x' . 
This pattern is the same as that in Ref. @, , m con- 
trast with x = x 1 found in Ref. Q. The density has a 
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FIG. 1: (Color online) Zero-temperature phase diagram for 
the triangular lattice. Second order phase transitions are de- 
noted by solid lines, whereas us and /j,sf are spinodal lines 
as a function of inverse coupling strength t = t/V. The thick 
lines dashed line at fi — 3 indicate first order transition be- 
tween SS A to SS B. Inset shows y/3X\/3 ordering of the solid 
and supersolid phases. 

discontinuous jump at Jl = 3, hence the SS A - SS B tran- 
sition is first order. Larger cluster size (9 sites) confirms 
this result. All these results (and more below) validate 
the SCMFT approach to the hard-core boson problem. 
The spinodal lines Jl (J) for the supersolid phases (not 
shown) have roughly a parabolic shape, closing at the 
critical endpoints (Jl = 3, t = 0) and (Jl = 3,t c = 0.216), 
the latter being the SS to SF transition. The maximum 
size of the metastable region, Jl = ±3.01, occurs halfway 
between t = and t c . 

The main properties of the supersolid phases are sum- 
marized as follows at the particle-hole symmetric point 
Jl = 3 (half-filling) . When t approaches 0, the super- 
solid state is in close proximity to the insulating states 
p — 2/3 (p — 1/3), therefore the jump in density dp 
is maximum in this region. The staggered density M s 
is also maximum there and vanishes continuously at the 
critical point t c = 0.216 after which only superfluidity 
survives. The superfluid density p s corresponds to the 
spin-stiffness in spin language. It measures the energy 
cost to introduce a twist 6 of the direction of spin between 
every pair of neighboring rows. We use its generalization 
to finite temperature following Ref. [l^- The SS A to 
SF transition is a continuous quantum phase transition 
with a kink in p s at the transition point. The value of p s 
that we find there (0.18) is within a few percent of the 
QMC results @. 

Finite temperature phase diagram at finite doping. In 
Fig. 2 we present the finite-temperature phase diagram 
along a vertical line t = 0. 1 of Fig. 1 . Because of particle- 
hole (Ising) symmetry it is a sufficient to show Jl > 3. 
Over a wide range of chemical potentials at high tem- 
peratures, a first order S to NF phase transition (dashed 
line) ends at a tricritical point d at about Jl = 3.70, where 
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FIG. 2: (Color online) (a) p-T Phase diagram corresponding 
to a vertical line t = 0.1 in Fig. 1. The inset shows the 
behavior near the particle-hole symmetric point p = 3. The 
two arrows indicate the region of metastability of supersolid 
phases between T CQ and T c . 
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FIG. 3: (Color online) (a) Density (p) hysteresis and super- 
fluid order parameter (*&) as a function of T along the p = 3 
vertical line of the phase diagram in Fig. 2. (b) Same quanti- 
ties but this time as a function of p along the horizontal line 
T — 0.06 in the phase diagram of Fig. 2. 



second order melting transition of the solid begins. Point 
e at the other end of the first order line marks the be- 
ginning of a very interesting region at large p. The first 
order transition bifurcates: to the right into a (BKT) 
transition separating SF and NF and, to the left, into a 
first order transition separating SF and S. Between point 
e and point f , we find the remarkable sequence of phases 
described in the introduction: As we raise the temper- 
ature at fixed p, one encounters SF, S then NF. The 
superfluid solidifies as we increase temperature because 
of an analog of the Pomeranchuk effect, the role of spin 
entropy being played by hard-core boson occupation of 
optical lattice sites. Solidification does not quench all 
the entropy. Let us come back to the BKT transition 
to the right of point e. One does expect the SF to NF 
transition to be of this nature 



14|. Clearly, SCMFT can- 




not accurately describe the topological BKT transition. 
Nevertheless, we take the jump in superfluid density p s 



FIG. 4: Density as a function of p for various values of T, 
each of which corresponds to a different horizontal cut on the 
phase diagram of Fig. 2. The inset shows the BKT transition 
in p s as a function of T at p = 3.2 and 6.2. The straight 
diagonal line is the BKT prediction. 



illustrated in the inset of Fig. 4, and the continuous van- 
ishing of the order parameter '5 = (S x ) = ((£> + /2) 
as very clear SCMFT signatures of the BKT transition. 

Supersolid phases appear near the symmetric point 
ju = 3. For pt < 3.38 the solid freezes into various super- 
solid phases with decreasing temperature. For example, 
at p = 3.2 the staggered density M s and the density p 
change continuously from S to the finite pt extension of 
SS A, but again there is a jump in the p s , as shown the 
inset of Fig. 4, so the transition is of the BKT type. 

The inset in Fig. 2 is a blow up of the region around 
the particle-hole symmetric point p = 3 where super- 
solid phases appear. At p = 3, raising T from zero, 
we notice that the ordering pattern of the solid changes 
from "ferrimagnetic" SS A (x, x, —2x') to "antiferromag- 
netic" SS C (x, —x, 0) at T Co = 0.053, indicated by point 
c in the inset. The SS A to SS C transition is first 
order, as can be seen from the hysteresis in the plot 
of density as a function of T in Fig. 3(a). The region 
of metastability associated with this transition is in the 
range f c (= 0.043) <f< f Co . The SS C phase contin- 
ues to higher temperature, up to point b. The SS C to 
solid transition point b (of BKT type) is indicated by 
the second arrow in Fig. 3(a). The area delimited by 
cO cl b contains the ferrimagnetic (x, — x', x") su- 
persolid phase that evolves from SS C with increasing pt 
for f Co <f < f b . The cO - cl line is second order. De- 
pendence on pt at fixed T = 0.06 for the superfluid order 
parameter and the density is shown in Fig. 3(b). The 
arrow to the right indicates the BKT transition from SS 
A to S: In Ref. 0, a BKT transition from SS A to S 
was also found with QMC at very similar temperatures. 
The arrow to the left marks the transition from SS C 
to SS A. The region b cl c2 delineates the solid or- 
der that evolves from 0) pattern. Outside this 
region, the solid phase has ordering (x,x,—x') and the 
transition between the two types of solids is second order. 

The plot of density as a function of pt in Fig. 4 confirms 
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FIG. 5: (a) The ground-state phase diagram for the t — t' — V 
model at p — 3 + . The solid line indicates second order tran- 
siton, whereas broken line indicates the first odrer transition. 
t' and V are measured with respect to t. The inset shows p 
vs t! hysteresis curve at V = 6. (b) M s , p — 1/2, and $ as a 
function of t' for V — 6. 

the order of the last two transitions we mentioned. The 
first kink in the T = 0.4 and that in the T = 0.1 curves 
are associated with, respectively, the second order solid 
to solid and SS C to SS A (at finite Jl where both phases 
are ferrimagnetic) . 

Finite t'. Finally, we investigate whether second- 
neighbor hopping t' , in the particle- hole symmetry case 
ju = 3, can induce the SS C phase at zero temperature, 
as proposed in Ref. [12]. A finite t' allows same sublat- 
tice hopping. In the presence of t', we choose clusters 
'1' and '2' shown in the inset of Fig.l, and connect them 
to each other through the perturbation t'. The effect of 
the other bonds that connect clusters T' and '2' are in- 
cluded in the self-consistent Eqs. (j4]). We checked that 
the ground state energy of this cluster is lower than that 
of cluster '3' (where all bonds reside on the cluster). 

The ground-state phase diagram for jus = 3 + is shown 
in Fig. E] (a). We note that for V = V/t > 3.0 a small 
value of the perturbation t' = t'/t drives SS A to SS 
C through a strong first order transition, as can be seen 
from the hysteresis exhibited in the inset of Fig. [5] (a). In 
part (b) of the same figure, we plot the staggered density 
M s , the superfluid order parameter 'F, and the average 
value of p— 1/2 as a function oft' at Jl = 3, corresponding 
to an horizontal cut at V = 6 in the phase diagram. We 
note that the finite value of p — 1/2 corresponds to SS 
A. With increasing t' the value of p — 1/2 jumps to zero, 
indicating SS C. Similar jumps can be seen in the other 
two curves. 

In summary, the strong geometric frustration present 
on the triangular lattice has striking consequences on the 
phase diagram of hard core bosons. First, as is well 
known, it allows the "ferrimagnetic" supersolids SS A 
and SS B phases to appear at T = 0. Second, the tri- 
angular lattice is associated with strong entropic effects 
at finite T that, we have shown, lead to a pronounced 
Pomcranchuk effect. We have also shown at the particle- 



hole symmetric point Jl — 3 that entropic effects at finite 
T, or finite t' at T = 0, lead to the appearance of the 
elusive "antiferromagnetic" SS C phase. Since the SS A 
and SS B supcrsolids break particle-hole symmetry, it is 
natural that increasing temperature restores a symmet- 
ric SS C state. In the case of t', it is a simple exercise 
to show that for same sublattice hopping, kinetic energy 
is minimized by (|0) + |1)) /y/2, i.e. the state in spin 
language. Finite t' thus also favours the restoration of 
the SS C (x,-x,0) state. The SS A to SS C transition 
is strongly first order under the influence of cither T or 
t' at Jl = 3 + . It is clearly not possible to see non-Landau 
Quantum Critical Point [HI with SCMFT, nevertheless 
it is likely that transitions that are strongly first order in 
SCMFT will not become continuous unless quantum fluc- 
tuations beyond the cluster size are singular enough to 
completely drive the transition. This is a delicate point 
that requires much more detailed studies guided by our 
results for phase boundaries. Our finite temperature re- 
sults are important for experimental studies of this very 
rich phase diagram with optical lattices or in solid state 
XXZ spin analogs. 
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